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Photographic  observations  of  the  Martian  satellites  were  made  at  the  opposition 
of  1967  with  the  Naval  Observatory’s  61-inch  astrometric  reflector.  A  small 
partially  transparent  metallic  film  filter  was  used  to  dimmish  the  light  from 
Mars  in  order  that  a  measurable  image  for  the  planetary  disk  as  well  as  for  the 
satellites  could  be  obtained.  The  plates  were  reduced  by  the  method  of  plate 
constants  using  positions  for  the  faint  background  stars  determined  from  astro - 
graphic  field  plates.  The  random  mean  error  of  these  observations  was  estimated 
to  be  not  greater  than  +0'T0. 

The  main  result  of  the  orbital  adjustment  is  a  +2°  correction  to  the  zero  of 
mean  longitude  for  Phobos.  This  confirms  the  findings  of  Wilkins  (1970)  and 
is  compatible  with  the  results  of  the  Mariner  9  observations.  The  scale  of  the 
orbits  of  both  satellites  gave  accordant  values  for  the  mass  of  Mars  and  the 
combined  value  of  30  99500  ±2800  (m.e.)  is  in  good  agreement  with  modern 
determinations. 

The  mean  error  for  Deimos  derived  from  the  residuals  after  solutionis  +0T11, 
which  agrees  well  with  the  observational  error  and  indicates  no  large  systematic 
error  in  either  the  theory  or  the  observations.  Por  Phobos,  however,  the  residual 
error,  ±0"19,  is  twice  the  expected  observational  error.  The  implications  of  this 
discrepancy  are  discussed. 


Introduction' 

In  liis  analysis  of  all  available  observa¬ 
tions  of  the  Martian  satellites  Wilkins 
(1964,  1967,  1968)  could  neither  confirm 
nor  reject  the  secular  accelerations  hi  the 
longitudes  of  Phobos  and  Deimos  found 
by  Sharpless  (1945)  20  years  earlier. 
Although  a  sizable  acceleration  term  was 
obtained  in  his  formal  solution,  Wilkins 
was  reluctant  to  accept  its  validity  because 
the  standard  error  of  the  solution  had 
increased  to  ±0"55  from  the  expected 
value  of  ±0"3  for  a  solution  using  only  the 
observations  from  a  given  opposition, 
telescope  and  observer  (visual).  The  ano¬ 
malous  increase  in  the  residuals  was  attri¬ 
buted  to  either  an  incomplete  theory  or 
systematic  '  observational  errors,  or  to 
both. 

1  Presented  at  I.A.TT.  Colloquium  No.  28, 
“Planetary  Satellites”,  Ithaca,  N.Y.  August, 
18-21,  1974. 
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While  a  reinvestigation  of  the  theory 
was  in  order,  and  has  in  fact  been  com¬ 
pleted  by  Sinclair  (1972),  there  are  at  least 
two  indications  that  a  large  portion  of  the 
anomalous  residuals  are  due  to  the 
observations.  First,  if  the  orbital  models 
are  too  simple,  it  is  expected  that  the 
effects  would  be  more  noticeable  for 
Phobos  than  for  Deimos.  This  is  not  the 
case;  the  standard  error  for  Deimos  is  a 
little  larger  than  that  for  Phobos.  Second, 
the  scales  of  the  orbits  of  the  two  satellites 
have  never  been  compatible  with  one  mass 
for  Mars  nor  with  modern  determinations 
of  the  mass. 

Late  in  1966  a  satellite- observing  pro¬ 
gram  was  reinitiated  at  the  Naval  Observ¬ 
atory  with  the  primary  aim  of  obtaining 
high-precision  observations  of  the  Martian 
moons.  To  date,  the  Martian  satellites 
have  been  observed  at  four  oppositions 
beginning  with  that  of  1967.  This  paper 
discusses  the  1967  observations. 
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Observational  Work 

The  superiority  of  the  photographic 
technique  for  the  Galilean  satellites  had 
been  demonstrated  early  in  the  century. 
Its  failure  for  the  moons  of  Mars  was  due 
to  the  proximity  of  the  satellites  to  the 
bright  primary  and  to  the  inferior  emul¬ 
sions  of  that  time.  Although  Sharpless 
made  a  considerable  effort  to  photograph 
the  moons  during  the  oppositions  of  1939, 
1941,  and  1943,  he  had  very  little  success. 
Good  photographs  of  the  satellites  were  not 
obtained  until  1954  by  Kuiper  (1961)  and 
others.  It  is  evident  that  a  long-focus 
instrument  is  needed  to  diminish  the 
intensity  of  the  halo  and  that  a  moderate 
/-ratio  is  necessary  to  insure  as  short 
exposures  as  possible  without  undue  loss 
of  astrometric  precision.  The  Naval 
Observatory’s  61 -inch  astrometric  reflector 
has  these  qualifications  in  addition  to 
having  astrometric  stability  and  was 
therefore  used  to  obtain  the  1967  series. 

The  early  success  of  the  photographic 
technique  in  satellite  astrometry  was  due 
largely  to  an  orbital  correction  technique 
in  which  only  intersatellite  measures  are 
used.  In  this  scheme  the  observed  and 
computed  positions  of  a  pair  of  satellites 
relative  to  each  other  are  used  to  correct 
the  parameters  of  both  orbits  simul¬ 
taneously.  Thus,  the  problem  of  making 
measurements  on  an  overexposed  primary 
is  eliminated.  This  method  however  is  not 
without  drawbacks.  The  most  obvious 
one  is  that  there  are  now  twice  as  many 
unknowns  in  each  conditional  equation  to 
determine  from  the  observations.  This 
alone  'will  degrade  the  precision  of  the 
corrections  but,  in  addition  to  this, 
statistical  correlations  between  the  correc¬ 
tions  to  the  parameters  of  both  satellites 
will  further  confound  the  results.  Another 
problem  is  encountered  when  the  observa¬ 
tions  for  one  satellite  are  more  accurate 
than  those  for  the  other  or  when  the  theory 
for  one  of  the  satellites  is  more  complete 
than  that  for  the  other.  This  is  the  case  for 
the  Martian  moons — the  observational 
and  theoretical  problems  are  greater  for 
Phobos  than  for  Deimos.  The  result  of  an 
orbital  adjustment  for  Deimos  from  obser¬ 


vations  relative  to  Phobos  would  not  lead 
to  an  improvement  in  its  orbit.  Further¬ 
more,  neither  the  ( 0  —  G)’s  nor  the  residuals 
after  solution  would  give  an  indication  as 
to  whether  the  problem  was  observational, 
theoretical,  or  for  that  matter,  which 
satellite  was  in  error. 

A  more  serious  complication  encountered 
in  the  use  of  intersatellite  positions  was 
pointed  out  by  H.  Struve  (Laves,  1938). 
Struve  showed  that  for  evenly  distributed 
observations  the  eccentricity  and  apse,  as 
determined  from  intersatellite  positions, 
have  only  1/30  the  weight  of  those  obtained 
from  planet-satellite  positions.  This  is 
important  since,  if  the  eccentricity  and 
apse  are  not  well  determined,  then  neither 
is  the  semimajor  axis ;  thus  the  mass  of  the 
primary  is  also  poorly  determined. 

It  is  clear  from  these  arguments  that 
observations  of  the  satellites  of  Mars 
should  be  measurements  relative  to  the 
primary.  In  order  to  obtain  a  measurable 
image  of  the  primary  a  technique  suggested 
to  the  writer  by  L.  W.  Frederick  (personal 
communication,  1966)  was  used:  A  small 
partially  transparent  metallic  film  filter 
“spot”  was  deposited  by  evaporation  onto 
the  center  of  a  GG14  filter.  The  metallic 
film,  winch  is  chromium  or  nichrome, 
transmits  neutrally,  is  microscopically 
thin,  and  does  not  affect  the  resolution.  Its 
optical  density  is  about  3.0,  which  equalizes 
the  intensity  in  the  image  of  the  primary 
to  that  in  the  satellites.  In  the  observing 
procedure  the  small  metallic  film  filter  is 
placed  over  the  primary.  [For  additional 
details  and  photographs  see  Pascu  (1967)]. 
The  plates  are  103aJ  in  combination  with 
the  above-mentioned  Schott  GG14  filter. 

For  the  1967  opposition,  69  positions 
were  obtained  for  Phobos  and  78  for 
Deimos,  on  six  nights  and  distributed  over 
three  weeks.  The  plates  were  reduced  by 
the  method  of  plate  constants  using  the 
faint  background  stars.  Positions  for  these 
faint  stars  were  obtained  from  astrographic 
field  plates. 

The  precision  of  the  observations  was 
estimated  by  two  methods:  First,  a 
comparison  of  independent  reductions  of 
the  direct  and  reverse  measurements — 
accounting  for  the  systematic  portion — 
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gave  a  mean  error  of  one  satellite  position 
with  respect  to  Mars  as  ±0''085.  This 
method  also  showed  that  the  position  of 
Mars  could  be  obtained  with  the  same 
precision  as  that  of  its  satellites.  The 
second  method  is  more  direct  and  involves 
the  comparison  of  ( 0  —  G)’s  for  two 
exposures  separated  by  only  a  few  minutes. 
It  is  assumed  here  that  the  variation  in 
two  consecutive  exposures  is  due  to 
observational  error  only.  In  this  way  the 
mean  error  of  one  observation  for  Phobos 
was  found  to  be  ±0''08  and  for  Deiruos 

±orio. 

Orbital  Adjustment 

In  the  classical  model,  introduced  by 
Struve  (1898),  the  satellites  move  in 
Keplerian  orbits  which  are  inclined  at  a 
constant  angle  to  the  Laplacian  or  fixed 
plane.  The  nodes  of  the  orbits  on  the 
Laplacian  plane  regress  at  a  constant  rate 
while  the  longitude  of  the  apse  increases  at 
approximately  the  same  rate.  The  orbital 
model  is  defined  by  the  eleven  parameters 
listed  in  Table  I  together  with  their 
initial  values.  These  values  are  taken  from 
Wilkins  (1964)  except  for  e,  K,  and  rrR 
which  are  from  Sinclair  (1972). 

In  polar  coordinates  the  conditional 


equations  for  the  differential  correction 
procedure  are 

where  zf<q  is  the  correction  to  the  ith 
element,  dpldei  and  cte/de;  are  the  partial 
derivatives  of  the  position  angle  p  and 
separation  s  with  respect  to  the  ith 
elem  ent,  and  n  is  the  number  of  parameters 
of  the  model.  The  right-hand  members, 
or  the  ( 0  —  G)’s  in  s  and  p  are  simply  the 
differences  between  the  observed  and 
computed  (predicted)  values  of  the  co- 
dinates.  Analytical  expressions  for  the 
partials  have  been  derived  previously 
(Pascu  1972). 

From  observations  obtained  in  such  a 
brief  interval  of  time  (three  weeks), 
corrections  to  the  motion  in  longitude  and 
the  motions  of  the  node  and  apse  cannot 
be  determined.  Furthermore,  the  small 
motion  of  the  node  during  these  observa¬ 
tions  indicates  that  large  correlations  will 
exist  between  the  corrections  to  J x  and  j 
and  the  corrections  to  0oand  N:.  In  order  to 
derive  as  much  information  as  possible 


TABLE  I 

Assumed  Initial  Values  foe  the  Orbital  Parameters  of  the  Satellites  of  Mars  (Equinox 

and  Equator  of  1950.0) 


Epoch  JED 

Phobos 

(2414600.5) 

Deimos 

(2414800.5) 

1.  a 

Apparent  semimajor  axis  at  unit  distance 

12T91 

32"36 

2.  e 

Eccentricity 

0.0184 

0.0020 

3.  n 

Daily  mean  motion  in  longitude 

1128?8443 

285?16192 

4.  Iq 

Mean  longitude  at  epoch 

227?1 

333?87 

5.  N± 

Longitude  of  node  of  fixed  plane  on  equator 

46?9 

46?40 

6.  Jx 

Inclination  of  fixed  plane  to  equator 

37?57 

36°64 

7-  j 

Inclination  of  orbital  plane  to  fixed  plane 

0?9 

1?80 

8.  e0 

Argument  of  node  of  orbital  plane  on  fixed  plane 

80° 

358?3 

9.  K 

Daily  motion  of  node  of  orbit  on  fixed  plane 

— 0?4354 

— 0?01813 

10.  7T0 

Longitude  of  pericenter  at  epoch 

211° 

300° 

11.  TTr 

Daily  motion  of  pericenter 

+0?4345 

+0?01814 
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TABLE  H 


Epoch.  JED 

Phobos 

2414600.5 

Deimos 

2414800.5 

a 

12f927±0r011 

32f341  4-  0''010  (m.e.) 

m-ia 

3097700  ±  7000 

3099900  ±  3000 

e 

0.0225  ±  0.0008 

0.0026  ±  0.0004 

AIq 

+2?10  +  0?1 

— 

ir0 

— 

341°  +  5° 

rms  ( 0  —  (7)before  solution 

±  Of  31 

±  ori4 

rms  ( 0  —  0)  after  solution 

±  ori9 

+  Of  11 

Estimated  obs.  mean  error 

+  oro9 

±oro9 

Combined  mass 

3099500  +  2800 

Fly-by  mass  (Mariner  IV)  (Null  1969) 

3098708  ±  9 

“  Mass  of  Mars  in  reciprocal  solar  units. 


from  the  observations  and  to  test  the 
reliability  of  the  results  a  number  of 
solutions  were  made.  These  trial  solutions 
began  with  corrections  to  a,  e,l0,  N lt  J lt  j, 
80,  and  7 r0  as  unknowns.  Unknowns  that 
were  indeterminate  or  highly  correlated 
were  then  dropped  in  order  that  a  conver¬ 
gent  solution  be  obtained. 

For  Phobos,  the  root  mean  square  (rms) 
(0  —  C )  before  solution  is  ±0"31.  The  first 
solution,  involving  corrections  to  the  eight 
parameters  mentioned  above,  diverged 
since  the  corrections  to  j  and  80  were 
indeterminate  and  the  correction  to  v0, 
while  significant  with  respect  to  its  error, 
was  too  large  to  be  handled  by  a  linear 
approximation.  In  order  to  obtain  a  better 
starting  value  for  7T0,  nine  additional 
observations  from  1969  and  1971  were 
added.  While  the  accuracy  of  these 
observations  is  not  as  good  as  those  of 
1967,  their  weight  in  the  solution  for  Att0 
is  greater.  Corrections  of  —9°  for  7r0  and 
—13°  for  0o  were  obtained  from  a  convergent 
solution  using  these  additional  observa¬ 
tions  and  including  corrections  to  all 
eleven  parameters  of  the  orbit.  A  conver¬ 
gent  solution  was  obtained  for  the  1967 
observations  -with  corrections-  to  a,  e,  l0, 
Nu  and  J x  using  the  corrected  values  for 
7 r0  and  60.  Significant  corrections  were 
obtained  for  a,  e,  and  l0  only  and  are  given 
in  Table  II. 

For  Deimos  the  situation  is  different. 


The  rms  ( 0  —  G)  before  solution  for  the 
1967  observations  is  only  ±0"14.  In  the 
orbital  adjustment  only  corrections  to 
a,  e,  and  tt0  were  determinate  and  the 
solution  with  these  three  parameters 
converged  rapidly,  giving  the  values  in 
Table  II. 

Conclusions 

There  are  a  number  of  conclusions  which 
can  be  drawn  from  the  results  presented  in 
Table  II  concerning  the  accuracy  of  these 
observations.  First,  the  good  agreement 
between  the  recomputed  rms  (O  —  C)  after 
solution  and  the  estimated  observational 
error  for  Deimos  indicates  that  there  is  no 
significant  systematic  error  in  the  mea¬ 
sured  position  of  the  primary  [that  is,  there 
is  no  center-of-figure  error  (Pascu,  1975)]. 
Second,  the  agreement  of  the  values  of  the 
mass  of  the  primary,  obtained  independ¬ 
ently  from  the  orbital  scales  of  each 
satellite,  and  the  fly-by  value  indicates 
that  there  are  no  significant  effects  due  to 
the  photographic  emulsion  or  to  the  halo 
of  the  planet. 

The  correction  of  +2°  to  the  longitude  of 
Phobos  confirms  the  value  found  by 
Wilkins  (1970)  for  1967  and  shows  that, 
despite  claims  to  the  contrary,  Earth- 
based  observations  can  and  have  in  fact 
detected  the  longitude  error  for  Phobos.  A 
longitude  error  of  this  magnitude  tends 
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to  support  the  secular  acceleration  theory 
for  Phobos  although  conclusive  evidence 
of  this  can  only  come  from  a  program  of 
concentrated  observation  over  a  period 
of  10  to  25  years. 

The  fact  that  the  residuals  after  solution 
(recomputed  ( 0  —  (7)’s)  were  comparable 
to  the  observational  error  for  Deimos 
suggests  that  Struve’s  theory  is  sufficient 
to  describe  its  motion  even  with  the 
considerable  increase  in  observational 
accuracy.  On  the  other  hand,  the  somewhat 
larger  than  expected  value  for  the  residuals 
after  solution  for  Phobos  means  either  that 
one  of  the  highly  correlated  parameters 
(for  the  1967  observations  only)  of  the 
theory  needs  a  substantial  adjustment  or 
that,  in  fact,  these  observations  are  precise 
enough  to  warrant  the  application  of  a 
more  sophisticated  theory,  such  as  that  of 
Sinclair  (1972). 
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